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Abstract 



The scattering of vortices at a critical value of the coupling constant in the Lagrangian can 
be approximated by a geodesic motion in the moduli space of classical static configurations of 



vortices. In this paper we give a scheme for generalising this idea to couplings that are near to 



the critical value. By perturbing a critically coupled field, we show that scattering of vortices 
at near-critical coupling can be approximated by motion in the original moduli space with a 
perturbed metric, and a potential. We apply this method to the scattering of two vortices, 
and compare our results to recent numerical simulations, and find good agreement where the 
scattering is not highly sensitive to radiation into other field modes. We also investigate the 
^ ■ possibility of bound stable orbits of two vortices in the quantum field theory. 
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1 Introduction 



Solitons in two and three space dimensions are well-known products of theories with spon- 
taneously broken symmetries. Examples are vortices in two dimensions, and cosmic strings, 
monopoles and skyrmions in three dimensions [0. Individual solitons are solutions of the clas- 
sical field equations, with a localised energy density and conserved topological charge. 

A difficult task is to determine the behaviour of soliton-soliton interactions, especially for 
theories where the solutions of the field equations cannot be written in closed form, but must 
be computed numerically. One approach that has attracted study is to put the field equations 
on a lattice, and to simulate the field theory on a computer. This approach has been fruitful in 
studying the interactions of vortices and strings , and monopoles and skyrmions to a lesser 
extent, but is however limited to describing the behaviour for definite values of parameters of 
the theory, and cannot provide a general description of soliton scattering. 

Another approach has been to examine the configuration space of the field theory (or rather, 
the "true" configuration space with gauge transformations quotiented out). The aim is to deter- 
mine a finite dimensional subset of this space on which the field equations can be approximated 
by a simpler particle-like dynamical system. Motion perpendicular to this space is ignored in 
the expectation that it doesn't have a significant effect on the dynamics for low velocities in 
the scattering. The scattering of monopoles || and vortices || has been modelled by adiabatic 
motion through the space of static configurations of the field theory at fixed topological charge, 
known as the moduli space. In the case of monopoles and vortices, current theory has only been 
able to model the interactions at a certain value of a dimensionless coupling constant in the 
Lagrangian, known as the critical coupling. This is known as the moduli space approximation, 
and was first proposed by Manton 0. 

This paper describes how the moduli space approximation can be extended to describe 
soliton scattering at near-critical coupling. Although our results are generalisable to other 
theories, we concentrate on the scattering of vortices in the Abelian Higgs model in (2 + 1) 
dimensions. This is a generic example of soliton scattering and describes the interactions of 
straight, parallel U(l) cosmic strings or magnetic flux tubes in a superconductor in the Landau- 
Ginzburg model. Sec. 2 outlines the theory of the moduli space approximation at critical 
coupling, and its extension to near-critical coupling. Sec. 3 tests this theory by applying it to 
the scattering of two vortices, and compares the results with numerical simulations of the full 
field theory. We conclude with a summary of our results. 

2 Extending the moduli space approximation 

We work with the Abelian Higgs model, with the vortices defined on R 2 , with a Lorentzian 

metric of signature (H ). [] The model has spontaneously broken U(l) gauge symmetry, 

and action 

S = J afxl—J^r + ^(KD^y + a(\<P\ 2 - c 2 ) 2 ] (1) 

1 The vortices could also be denned on a general orientable Riemann surface such as S 2 , S 1 x S 1 or hyperbolic 
space H 2 . We take K 2 here for simplicity. 
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where is complex scalar valued. D^cp = (c^ — ieA^)(p where is the real valued gauge 
potential, and e the gauge coupling. The action (|]) can be rescaled using the dimensionless 
quantities 



T /" _ ; 

JU v 

ec 

(f) = Cq4> 

= c A 



(2) 



S = 



to read 



s = J d 3 x[-\f^r + + (\ + A)(|0| 2 - 1) 2 ] (3) 

where | + A = ^ (so A > — ~) and D^cfi = (<9 M — iA^)(j). In these units the mass of the gauge 
particle m p hoton = 1, and the mass of the scalar field is m Higgs = \/T+ 8A. We work with this 
action, dropping the tildas. The action can be split into kinetic and potential pieces, which 
after choosing the gauge Aq = can be written 

S = J dx (T - V) (4) 

T = J d 2 x^(4>4>* + AiAi) 

V = /^[^/i 2 2 + ^A0(A0r + (^+A)(|0| 2 -l) 2 ] 

(the equation of motion corresponding to A must then be imposed as a constraint; Gauss' 
Law). T also provides a natural metric on configuration space. Any finite energy configuration 
has a magnetic flux 

" d 2 xf 12 = 2ttN (5) 



where N is an integer, and this magnetic flux is localised in the neighbourhoods of the N zeros 
of 4>. This topological identity can be used to write the potential energy as a sum of squares, 
plus a remainder term. Introducing the complex coordinates z = x\ + 1x2, z = X\ — ix 2 , 

V = vriV + / d -^[\{h, + ^(I0| 2 - I)) 2 + \\D- Z <P? + A(|0| 2 - I) 2 ] (6) 

At A = the last, remainder, term in the above equation vanishes. The potential energy V can 
then saturate its topological lower bound Ntt when the Bogomolny equations M 

D s <j> = (7) 

/l2 + ^(|0| 2 -l) = 

are satisfied (A = is the critical coupling for this theory). The first equation of (0) can be 
solved to give 

A- z = -id- z \og<P (8) 
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so that the gauge field is dependent on the Higgs field, and can be eliminated as a variable 
in the equations (0). When T = and V = Ntt, these equations give the static solutions of 
the field equations (cj)^(Q, z, z), A\ \Q, z, zj). Q is a set of collective coordinates representing 
the position of the field in the moduli space. This is because at critical coupling there are no 
forces between static vortices. It has been shown J7J that on any orientable Riemann surface, Q 
can be specified by the unlabelled set of the zeros of (f>(°\ {z 1 . . . z N }, which may be identified 
with the positions of the vortices. An important aspect of soliton behaviour is that solitons are 
indistinguishable particles at the classical level, and so the coordinates Q are invariant under 
the action of the permutation group which exchanges vortex positions. The moduli space 
of iV vortices defined on an orientable Riemann surface M is then Airv = 

The essence of the moduli space approximation at critical coupling is to project the dynamic 
field evolution onto the moduli space. The projection <p — > <p^\Q) is made by fixing the zeros 
of the Higgs field, and the dynamics of the field are approximated by allowing the coordinates 
Q to vary in time. The time derivative of the Higgs field is approximated by 

, (o) = lim ^m+5t))-^m)) 

y st-,o 8t 

and similarily for the gauge field A^. Samols |4[] has shown that, using the approximation (|9]), 
the kinetic energy of the fields can be written as 

N 1 dzdz- 
r = V -n 9l AQY-^-^r (10) 
.f- 2 yi3WJ dt dt K J 

9ij = 6 ij+ 2- 

where bj is the linear coefficient in an expansion of f(z, z) = log l^*- -*! 2 about the vortex at zf 

f(z,z) = log\z-z j \ 2 + a J {Q) + h J {Q){z-z J ) + h j {Q){z-z J )+ (11) 
Ci(Q)(z - zj) 2 + Cj(Q)(z - zj) 2 + dj(Q)(z - Zj)(z — Zj) + . . . 

When there are no massless fields present, as in this theory, motion perpendicular to the moduli 
space is suppressed by the positive quadratic nature of the potential in the directions orthogonal 
to .Mat. The normal frequencies are bounded below by the photon mass and Higgs mass. Whilst 
no rigorous field theory analysis has been done, finite dimensional theory suggests that in a 
scattering process where the timescale of the process along the moduli space is t± and the 
timescale perpendicular to the moduli space is t 2? the fraction of energy transferred to modes 
perpendicular to the moduli space is 0(exp —ti/t 2 ). For vortex scattering with typical velocity 
v , this is 0(exp —-) and so is suppressed at small v. This approximation has been tested by 
computer simulations, and found to work well up to velocities v ~ 0.4 [|J. 

Expressing the motion of solitons in terms of a drift of a projection of the field onto the 
moduli space of static configurations suggests a similar approach when A ^ 0, A small. When 
A 7^ 0, the vortices exert a static force on each other, and the space of static solutions collapses 
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to isolated points in configuration space when spatial translations and rotations are quotiented 
out (for example, all N vortices coalesce into a stable winding number N vortex in the attrac- 
tive case). One possible way to model the interactions of weakly attracting/repelling vortices 
would be to project the motion onto the set of curves of steepest descent of the potential, 
connecting unstable to stable static configurations. This technique originates in molecular re- 
action dynamics, and might be applied to the scattering of skyrmions ||. The true evolution is 
expected to follow a path in some neighbourhood of these gradient curves, for sufficiently small 
kinetic energies. However, the calculation of these gradient curves is a difficult constrained 
optimisation problem. In this paper, we argue that the dynamics is equally well approximated 
by projection onto a simpler, nearby manifold: the moduli space at critical coupling. 

We conjecture that as the potential energy of a particular configuration V[<f>, A{, A] changes 
smoothly with A according to equation (^), the set of steepest descent curves lies in the neigh- 
bourhood of the moduli space of static configurations (4>(°\ A\ ^) at critical coupling (for suf- 
ficiently small A), and indeed is diffeomorphic to it. Although these latter configurations are 
now no longer static solutions of the field equations or of minimal potential energy, for A small 
enough the evolution of the fields will remain in their neighbourhood. 

We propose an approximate description of the scattering by projecting the field at any given 
time to A = as well as T = 0, keeping the zeros of the Higgs field fixed. The scattering is then 
modelled by collective coordinate motion in the original moduli space at critical coupling. The 
induced metric on this space is a perturbed version of the metric at critical coupling, however 
we give arguments below that it is not necessary to accurately determine this perturbation for 
all vortex positions. There is also an induced potential on the moduli space. We may write 

(f) = (J>( )(Q(t))exp8h(z,z,t) (12) 
A = 4°\Q(t)) + 8Bi(z,z,t) 

where h and B^ are 0(1), but 5 is small. Taking e to be a measure of the (small) kinetic 
energy of the fields, we assume that for A and e sufficiently small, S(X, e) is 0(A) small. There is 
actually a technical complication with ( |T2"D as h and Bi are not independent, but are constrained 
by Gauss' Law. Furthermore, because of the residual time-independent gauge symmetry, the 
solution may evolve to be close to a gauge transformed multivortex. We will not need to take 
account of this as a knowledge of h and Bi are not necessary for our approximation. The point 
of equation ([12]) is to estimate the nearness of the fields to the moduli space at a particular 
time, and their kinetic energy. Substituting (12) into (|6|) we find that the potential 



V = Ntc + XJ d 2 x{\^ 0) {Q)\ 2 -l) 2 + O(X 2 ,X5,5 2 ) (13) 
= A% + XV (Q) + 0(A 2 ) 

is independent to first order in A of the deviation of the field from M.^. This is expected, because 
M.n was at a minimum of the potential energy for A = 0. Choosing a static configuration 
(f>(°'(Q) at A = and "turning on" A results in the minimum of V for fixed Q being at an 0(A) 
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distance from cp^(Q), but the difference between the new value of the potential of 4> iff) given 
by fll3|) and the new minimum at fixed Q is 0(A 2 ). 

The situation with the kinetic energy is not so simple; it varies to first order in 5, as can 
be seen by direct substitution of ([Tj|) into To determine the first order variation of the 
kinetic energy, we would in principle have to solve for h and Bi by minimizing the potential 
energy at fixed Q. However, we now give arguments that in a low velocity scattering this is not 
necessary It is not our intention to give a rigorous proof of our intended approximation, but 
rather to base our argument on the physical properties of vortex scattering. 

If we consider the scattering of near-critical vortices with typical low velocity v at a particular 
time, then the rate of angular deviation from the critical coupling motion that can be attributed 
to the potential is 0(X/v). The potential is only appreciably different from zero over the 
timescale of the scattering t* = 0(1 /v) so the total angular deviation due to the potential 
can be estimated as 0(X/v 2 ). The rate of angular deviation from critical coupling motion 
that can be attributed to the perturbation in the metric is 0(Xv 2 /v). This acts for all times, 
including before and after the scattering. But if we can calculate the effect of the perturbed 
metric for times \t\ > t* when the vortices are well separated, and make some approximation 
to it for the times when the vortices are close then the effect of our lack of knowledge of the 
metric perturbation when the vortices are close is a total angular error of 0(A). Thus for slow 
motion the effect of the potential dominates over the perturbation in the metric, provided the 
perturbation is correct asymptotically. 

To determine how the metric changes with A when the vortices are well separated, let us 
first consider a static N = 1 solution (<f>^\Aj) at A ^ 0. This has potential energy which 
from ( |I~3"D can be written as 

V = 7r + A|rf 2 x(|0 (o) | 2 -l) 2 + O(A 2 ) (14) 
= 7t + AVi + 0(A 2 ) 

where (f)^ is the Higgs field of a single A = vortex. V\ is a constant, whose numerical value 
is given in Sec. 3. The theory is Lorentz invariant, so we may obtain the energy of a moving 
vortex by a Lorentz boost, which in the low velocity limit becomes a Galilean boost. The 
energy is then 

E = (tt + XV 1 + 0(X 2 ))(1 + \ 2 + 0(v 4 )) (15) 
so we may identify the potential energy of a single vortex as its inertial mass 

m = 7r + XV\ + 0(A 2 ) (16) 

Consider now a configuration of slowly moving, well-separated vortices at near-critical coupling, 
such that the significant contribution to the kinetic energy comes from the motion of the vortices 
(rather than radiative modes). The Higgs and gauge fields are massive, and exponentially 
approach the vacuum away from the vortices. Thus the significant part of the kinetic energy 
comes from a region close to the zero of the Higgs field. We assert that in this region, the 
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configuration and its time derivative is well approximated by the boosted N = 1 solution 
(cj)^\Aj) with the appropriate velocity, centred at the zero of the Higgs field. The kinetic 
energy will then approach the kinetic energy of N moving particles of inertial mass m 

T~-mV- (17) 

2 f-f dt dt y 1 

This models the asymptotics of vortex dynamics. We therefore choose as our approximation to 
the kinetic energy 

N 1 dz 1 dz j 

T(Q,Q) = E 2 {7r + XVl)9 ^ Q) lF^ (18) 

where the metric gij is defined in ([Hj). As a consequence of the exponential decay of static 
solutions, the metric quickly becomes flat as the vortex separation increases, and T(Q, Q) 
quickly approaches (pp. 

In summary, the evolution of a configuration of vortices in a scattering problem is expected 
to be well approximated by a an evolution of particles (the zeros of the Higgs field) according 
to the action 

S = J dt(T(Q,Q)-XV(Q)) (19) 

where T(Q,Q) is defined in (|TSD and V(Q) is defined in (p~3|) . 

In a recent paper ||, Stuart has established some rigorous results on vortex dynamics. 
An existence theorem is proved for an initial value problem in configuration space, where the 
initial data is 0(|A|) close to the moduli space, the initial field derivatives are 0(|A|2), and 
0(| A|) close to the translational zero modes (at critical coupling) on the moduli space. For a 
subsequent time T* = |A| - i the evolution continues to satisfy the same conditions as the initial 
data. The point on the moduli space of closest approach to the field configuration is 0(|A|a) 
close to a Lagrangian motion given by the metric on the moduli space at critical coupling flTO|) 
and the potential V(Q) from (|T3|) . This applies for all vortex configurations, not just for well 
separated vortices. However, the metric used is incorrect asymptotically, as the inertial mass 
of the vortices changes with A. This contributes to the 0(|A|a) drift from the Hamiltonian 
evolution on the moduli space (this can be seen immediately by considering the motion of a 
single boosted near-critical vortex in the context of this problem). We anticipate that Stuart's 



results would still hold if our metric (|18|) was used together with the potential to give the 
Lagrangian evolution on the moduli space. Furthermore, in a scattering problem where the 
vortices are well separated except for a time of 0(T*), we conjecture that Stuart's results may 
be improved to hold for all times with the use of our asymptotic metric. 

We now seek to test the validity of our approximation against "observation": numerical 
lattice simulations of the full field theory. 
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3 The metric and potential of two vortices 



For two vortices, the geometry of the moduli space is very simple. Using the notation of Samols, 
we decompose the moduli space as the product 

M 2 = C x M 2 ° (20) 

where, if vortices with masses m lie at z\ and z 2 then the centre of mass coordinate Z = 
\{z\ + z®) ^ can be quotiented out, leaving the system with a reduced mass of m/2 and 
relative coordinate ( = \{z\ — z 2 ) G M 2 °. The vortices lie at ±£ in these coordinates. We will 
find it convenient to use the polar coordinates 

C = ae id (21) 

where, due to the symmetry of the scattering under interchange of two vortices, the range of 
6 is [0, 7r). Hermiticity of the metric on M 2 ° (in addition to rotational symmetry) follows from 
the reality of the kinetic energy, and requires it to have the form 

ds 2 = F 2 (a)(da 2 + a 2 d6 2 ) (22) 

F 2 (a) can be computed from the solution 0(°)(<r) of the Bogomolny equations with zeros of 
on the real axis at £ = ±cr by the formula 

F*(a) = l + l±(ab) (23) 
a da 

where b is the coefficient in the expansion of log |0^°^(cr)| 2 around the vortex at £ = a (see (TO))- 
As explained by Samols, the two dimensional space M® has the form of a rounded cone when 
embedded in Euclidean R 3 . 

The potential on the moduli space away from critical coupling is clearly a function of a 
only. Recalling (fl3|) it is calculated by 

V{*)= [ M \cr)\ 2 -l) 2 d 2 x (24) 



R, 

The interactions between vortices are smooth when expressed in terms of a good global co- 
ordinate. Due to the symmetry of the field under the interchange of two vortices, ( is not 
well-defined as £ — > 0, but uj = C, 2 G (D is. For small a the leading order behaviour of V(a) 
is determined by the requirement that exists at to = 0; in fact by rotational symmetry, 

£ = at u = 0. If we define V 2 = V(0), V(u) = V 2 + 0(cu 2 ) and 

V{a) = V 2 + aa 4 + 0(a 6 ) for a < 1 (25) 

By numerical computation, V 2 ~ 13.63 and a ~ 0.36. For large vortex separation, the contri- 
bution to the integral (p4|) from the neighbourhood of one vortex is only slightly perturbed by 
the presence of the other. For the Higgs field (p^ of a single vortex at £ = a 

< 1 - l^l 2 < M(5)e' {1 ~ 5)a for any 5 > 0, \z - a\ > a (26) 
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Figure 1: A plot of V(a) for two vortices where the inter- vortex separation is 2a. 



Figure 2: A comparison of the stationary vortex interaction potential calculated by Jacobs and 
Rebbi for A = 0.08625 (upper curve) and A = —0.06375 (lower curve), scaled to our units, to 
the potential V(a) — V(oo) (dashed line) 



which implies 

\<P i0 \*)\ 2 = \^W°l\ 2 + 0(e- 2 ^) (27) 

so 

2Vi < V(a) < 2V 1 + 0(e- {1 - S)a ) (28) 

Numerically, 2V X ~ 10.44 < V 2 . Q 

The numerical computation of the metric and potential had two stages. Firstly, the field 
(a) was calculated from the Bogomolny equations using an initial approximation which was 
a superposition of two approximate 1-vortex solutions. This was relaxed using a multigrid 
algorithm on a Silicon Graphics Indigo workstation. Seven grids were used, the densest of 
which was 256 x 384. Convergence was typically quite rapid, requiring only a few sweeps over 
the range of grids. The derivative of In |(/>^(<j)| 2 (for the metric) at the vortex position was 
calculated by fitting a cubic spline curve along the real axis joining the two vortices. The 
integral V(a) given by equation (|24]) was performed by fitting bicubic surface splines to the 
lattice solution. This was done for values of a from to 8 at intervals of 0.02. Finally, a 
suitably smoothed spline curve was fitted to the functions F 2 (a) and V(a). A graph of the 
potential V is shown in Fig. [I]. Jacobs and Rebbi [II] have calculated the interaction energy of 
two stationary vortices at coupling constants (in our units) of A = 0.08625 and A = —0.06375. 
They use a functional ansatz for the Higgs field which is a weighted superposition of the Higgs 
field for two isolated vortices at ±er and a 2-vortex at the origin, plus a correction term which 
is expanded as a power series in \z\ and cos(arg z). The free parameters are the coefficients of 
the power series which is truncated at sixth order. They then optimise the potential energy 
numerically as a function of those parameters. A scaled comparison of their results with our 
potential is shown in Fig. |2]. 



4 Scattering behaviour of two vortices 

The scattering of two vortices at critical coupling has for some time been studied numerically 
by putting the field theory on a lattice 0. Recently, Myers, Rebbi and Strilka [T2| simulated 
vortex scattering on a lattice at couplings A = — ^,0,| (in their units, this corresponded to 
A = 1,2,4) [0. Here, we compare our results at the values of the coupling of A = — ^ and 
A = | with theirs. For a comparison of A = scattering with numerical simulation, see ref . |4j . 

2 This argument can be generalised to N vortices. If {z\ . . . zjv} are the vortex positions, then NV\ < V < 
NVi + 0{e-V- s >) where here a = min i¥:j (\zi - z 3 )\. NVi <V N = V(0) 
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Figure 3: Some sample trajectories in the centre of mass frame for vortices with A = | and 
initial velocity v = 0.3. The lowest impact parameter graph is at d = 0.01. The trajectory of 
only one vortex is shown. 



Figure 4: A comparison of scattering angles calculated by Myers et al. to values calculated 
using the moduli space approximation when A = |. The impact velocities shown are v — 0.1 
(+), v = 0.2 (o), v = 0.3 (x), v = 0.4 (*) 



We turn first to A = | (type-II) vortices. The vortices were initially placed at large x- 
separations on paths with velocities (±t> , 0) and an impact parameter of 2d. The equations 
of motion derived from our model Lagrangian flT9"|) were integrated numerically using a fourth 



order Runge-Kutta method with variable step size. In Fig. y we show some trajectories for 
initial vortex velocity of v = 0.3. 

Scattering angles for a range of impact parameters and velocities are shown in Table 1. 
As the interaction is repulsive, there is a critical velocity below which a head on collision will 
scatter the vortices elastically back through 180°, and above which the vortices will scatter 
through 90°. This is calculable from the conservation of energy in our approximation, giving 



J crit 



'WHM.MVJ (29) 

27T 



so for A = |, v cr u — 0.25. For initial velocities below v cr u the scattering angle increases smoothly 
with decreasing impact parameter. For velocities above, but close to v cr i t the scattering angle 
has a turning point at very small impact parameter, and scattering in the backwards direction is 
still possible. As the velocity increases further, this turning point disappears, and all scattering 
is in the forward direction. A comparison with the results of Myers et al. is shown in Fig. |j. 
The agreement is very good at low velocities, but less good as the velocity increases. This is 
due to the non-relativistic nature of the moduli space approximation, and radiation into other 
field modes. This loss of quantitative accuracy is observed to a similar degree in the moduli 
space approximation for critical vortices. 

Turning now to A = — ^ vortices (type-I), the attraction due to the potential is longer 
range than the repulsion due to the metric (the mass of the photon is larger than the mass 
of the Higgs particle). Scattering trajectories for initial velocity v = 0.3 are shown in Fig. |5|, 
and Table 2 gives scattering angles for a range of initial velocities and impact parameters. The 
scattering behaviour is now qualitatively the same at all velocities. For large impact parameter, 
the scattering angle is small and negative. As the impact parameter decreases, the scattering 
angle decreases to a minimum, and then as the repulsion of the vortex cores overcomes the 
attraction of the potential, increases back through zero and approaches 90° scattering at zero 
impact parameter. Scattering at low impact parameter is insensitive to the vortex velocity. For 
initial velocity v = 0.3, and impact parameter d ~ 2.8 the repulsion of the vortex core balances 
the attraction of the potential, giving zero deflection. 
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Table 1: The scattering angle 9a e fi for A = | is plotted versus impact parameter and initial 
velocity v. 



impact parameter, d 



V = 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


0.05 


160.1 


140.3 


120.6 


101.2 


82.3 


64.0 


46.9 


31.7 


19.3 


10.5 


0.1 


151.7 


124.7 


99.5 


76.6 


56.1 


38.5 


24.2 


13.7 


7.0 


3.3 


0.15 


140.3 


106.9 


80.0 


57.7 


39.3 


24.7 


14.1 


7.2 


3.5 


1.5 


0.20 


119.1 


87.2 


63.5 


44.2 


28.7 


17.1 


9.2 


4.5 


2.1 


0.9 


0.25 


90.4 


70.8 


51.8 


35.4 


22.2 


12.6 


6.6 


3.1 


1.4 


0.6 


0.30 


77.3 


61.0 


44.3 


29.6 


18.0 


9.9 


5.0 


2.3 


1.0 


0.4 


0.35 


72.6 


55.7 


39.5 


25.7 


15.2 


8.1 


4.0 


1.8 


0.8 


0.3 


0.40 


70.5 


52.6 


36.4 


23.1 


13.2 


6.9 


3.3 


1.5 


0.6 


0.3 


0.45 


69.4 


50.7 


34.3 


21.2 


11.8 


6.0 


2.9 


1.3 


0.5 


0.2 


0.50 


68.7 


49.5 


32.9 


19.8 


10.8 


5.4 


2.5 


1.1 


0.5 


0.2 



Figure 5: Some sample trajectories in the centre of mass frame for vortices with A = — and 
initial velocity v = 0.3. The lowest impact parameter graph is at d = 0.01. The trajectory of 
only one vortex is shown. 



Table 2: The scattering angle Qdefi for A = is plotted versus impact parameter and initial 
velocity v. 













impact parameter, d 








V = 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


0.05 


66.6 


43.4 


19.1 


-7.0 


-36.1 


-70.1 


-113.9 


-182.6 


-72.6 


-12.8 


0.1 


63.7 


37.1 


8.8 


-22.9 


-60.7 


-109.4 


-61.0 


-16.9 


-6.2 


-2.5 


0.15 


63.2 


36.4 


8.5 


-20.2 


-42.2 


-32.7 


-14.2 


-5.8 


-2.4 


-1.0 


0.20 


63.6 


37.5 


11.9 


-10.4 


-19.5 


-13.4 


-6.5 


-2.9 


-1.3 


-0.5 


0.25 


64.0 


39.0 


15.4 


-2.7 
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Figure 6: A comparison of scattering angles calculated by Myers et al. to values calculated 
using the moduli space approximation when A = — A. The impact velocities shown are v — 0.1 
(+) (note there are few points here because of vortex capture) and v = 0.2 (o). 



Figure 7: As above, with v = 0.3 (x), v — 0.4 (*). 



Interesting behaviour is observed at low velocity. The largest (negative) scattering angle 
diverges as u — > 0; the vortices can orbit around each other, forming a "glory" (it has been 
shown [13 1 that this type of behaviour does not occur in the scattering of photons off cosmic 
strings, where the photons propagate in the background gravitational metric of the string). 
Unfortunately, as the scattering angle is very sensitive to changes in the impact parameter 
or velocity in this region, it is also very sensitive to radiation into other field modes, and 
so our approximation will only give a qualitative description of the scattering. Myers et al. 
observed that for low impact parameter and velocity, this radiation means the vortices can no 
longer escape their attraction, and scatter repeatedly through 90° with the maximum separation 
decreasing with each cycle; presumably the vortices would eventually settle down to a stable 
winding number two vortex. We cannot observe this end result in our approximation, as the 
conservation of energy and angular momentum forbids it, but we do observe bound orbits in 
the two vortex system. We turn to the question of their stability in the next section. Figs. ^| 
and [7] give a comparison our results to Myers et al. Agreement is best for mid-range velocities 
v = 0.2, 0.3, 0.4, but we do not expect our results to be quantitatively accurate for v — 0.1. 

A summary of the scattering angles versus impact parameters is shown in Fig. |8|, where the 
critical coupling result is also marked. 



5 Stability of bound orbits 

In this section we would like to consider the long term behaviour of a bound system of two 
vortices, from a classical and quantum mechanical point of view in the moduli space approxi- 
mation. The radial equation of motion for the two vortex system, expressed in the coordinates 
(a, 6) used in the previous section, can be rewritten using the conserved energy E = T + V 
(for convenience we normalise V here by subtracting V(oo) from it) and angular momentum 
h = a 2 F 2 6 as 

2rtF 2 & 2 = E-XV(a)-^^ = E-w(a,h) (30) 

Thus the problem is a slightly modified version of a particle in a potential well. For h ^ 0, the 
metric radial function F works effectively to reduce the angular momentum at smaller a and 
therefore increases the strength of the centrifugal barrier as a — ► 0. 



Figure 8: Summary of the scattering angles computed for A = (dashed line), A = | (above 
dashed line) and A = — A (below dashed line). 
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Figure 9: The effective potential w(a, h) for h £ [0, 0.1] at intervals of 0.01. 



When A < and h = 0, scattering with E > produces the familiar 90° behaviour. Motion 
with XV (0) < E < produces a "breather" state, with the vortices separating to a maximum 
distance, then returning and scattering through 90° in each cycle. For h ^ 0, analysis of the 
function w(a) shows that for h 2 < 1.14|A| , w(a,h) has two turning points, the minimum at 
radius a < 2.70, the maximum at radius a > 2.70. This implies the existence of trapped orbits. 
For h 2 > 1.14|A|, no orbits are possible. Fig. || illustrates the effective potential w(a, h) for a 
range of values of h. 

The minimum of w(a) corresponds to a stable circular orbit in our approximation. In the 
classical field theory, energy will be radiated from the system, and the vortices will spiral into a 
stable winding number two vortex. However, this decay might not occur in the corresponding 
quantum field theory. The available energy of the system is proportional to A, whereas the 
mass of the photon (smaller than the Higgs mass for our range of parameters) is proportional 
to the gauge coupling e, scaled to e = 1 here. Hence if A is sufficiently small, there will not be 
enough available energy for photons to be radiated. 

We do not attempt to analyse the quantum field theory here, but an alternative method 
suggests itself: treating the zeros of the Higgs field as quantum mechanical particles of mass m, 
and using Schrodinger's equation to give an approximate quantisation. For the wavefunction 
if), the time-independent Schrodinger equation is 



We separate of variables, by setting 

if; = e^R(a) (32) 

where I = 2nh , n G Z as 6 £ [0, ir). We can then reduce the resulting equation for R(a) to a 
standard form by the substitution R = x a ~* giving 

X "+ 2 -^F 2 (E-U)x = (33) 

where 

The wavefunction decays exponentially in regions where E < U(a, n) (if the vortex mass is large, 
as is typically the case for an artifact of symmetry breaking, the spectrum of allowed bound 
energy states is dense, and the system behaves semiclassically. The wavefunction is confined to 
the close neighbourhood of regions where E > U{a,n). However, we need not assume that the 
vortex mass is large here). As we are in two spatial dimensions, there is always at least one 
energy level with E < (contrary to the result for three dimensions) |14j]. For bound states 
with E > there is a small probability of tunnelling through the potential barrier, allowing 
the vortices to escape to infinity. States with E < are confined to a neighbourhood of the 
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classically accessible region. The maximum energy available for decay of a bound configuration 
with E < is \\\V dept h, where V depth = V(0) - V(oo). m Higgs = 1 + 0(A) is slightly smaller 
than m p hoton = 1, but this only changes the stability of an orbit to second order in A. Hence 
all configurations with E < are quantum mechanically stable if 

|A| < — ^ — ~ 0.31 (35) 
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